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In this paper, robustness of planar trivalent random graphs to targeted attacks of highest con- 
nected nodes is investigated using numerical simulations. It is shown that these graphs are relatively 

qq ■ robust. The nonrandom node removal process of targeted attacks is also investigated as a special 

case of non-uniform site percolation. Critical exponents are calculated by measuring various prop- 
erties of the distribution of percolation clusters. They are found to be roughly compatible with 

£NJ ■ critical exponents of uniform percolation on these graphs. 
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I. INTRODUCTION 

o , 

Robustness or fragility of a graph characterize its behavior when systematic or random deletion of a fraction of 
nodes is performed. These questions are of practical importance for real complex networks (World Wide Web, social 
networks, cells, etc) and, recently, there has been much interest in investigating targeted and random attacks for two 
families of random complex networks 0, 0, 0] , namely exponential and scale-free networks 0| • Resilience to random 
deletion of nodes can, for example, be seen as tolerance against error, i.e. random failure of a fraction of nodes. 
Targeted deletion of high degree nodes can simulate attacks of hackers on computer networks. However, besides these 
practical issues, investigating random and targeted attacks on abstract graphs is of great theoretical interest in its 
own right. It has to be noted that the random attack problem is equivalent to uniform site percolation 

In this paper, we study a family of purely abstract random graphs extensively used in the past decades as discrete 
models for euclidean quantum gravity (>\. namely random <£> 3 planar graphs. These planar graphs are made of trivalent 
vertices and look locally like the regular honeycomb lattice. But the faces are not necessarily hexagonal, so that long 
distance properties are very different from those of honeycomb lattices. Note that site percolation on these graphs 
has already been investigated @,@|- By means of Monte-Carlo simulations, we study resilience of these graphs against 
targeted attacks defined as a systematic face removal process, i.e. removal of all links and vertices that belong to 
the edges of the targeted face. More precisely, we investigate connected clusters distributions when all faces whose 
' number of edges is larger than a cut-off k max are removed. 

A key feature in the behavior of graphs against targeted attacks is the degree distribution. However, vertex 
CN ' coordination numbers on complex networks are uncorrelated random variables so that the degree distribution is 
sufficient to characterize these random graphs. On the contrary, planarity constraints induce long range correlations 
■ for random $ 3 planar graphs so that the degree distribution is not enough to characterize these graphs Q. This 
characteristic is also shared by Voronoi'/Delaunay graphs. However, correlations in the latter case decrease much 
faster than for random $ 3 planar graphs. This fundamental difference between complex networks and random $ 3 
planar graphs induces radically different behaviors under uniform percolations problems. The main purpose of this 
work is then to determine whether it also induces different behaviors against non-uniform percolation like targeted 
attacks, and, if so, to quantify these differences. 

It has to be noted that, even if random <f> 3 planar graphs have been used in Euclidean quantum gravity, the 
problem studied here is not relevant to these issues. Indeed, there is no back-reaction of any degree of freedom on 
the connectivity of the graphs. In other words, the disorder due to random degree distribution is quenched and 
we study the average behavior against targeted attacks of a graph picked at random in the ensemble of random $ 3 
planar graphs. However, this work could be the starting point of the study of a modified model of quantum gravity 
characterized by matter fields (Ising spin, for instance) coupled to the ensemble of graphs obtained from the random 
<I> 3 planar graphs after a systematic node delation procedure. 

Besides the size of the largest connected component S, we study several properties of connected cluster of graphs 
for various values of k max , in particular, the density of the second largest cluster S2 and the average size of finite 
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Figure 1: A $ 3 



planar graph 



clusters S. Considering k max as a continuous parameter, it is possible to extract a critical parameter k c similar to a 
percolation threshold, and to compute critical exponents (3, 7 and v associated to the second largest cluster, average 
cluster size and correlation length. These exponents are compared with those of uniform percolation problem on 
random <£> 3 planar graphs. 

This paper is organized as follows: in section UU we introduce the model of planar $ 3 random graphs; in section Ullj 
targeted attacks are defined and results of numerical simulations are reported; the problem is described as non-uniform 
percolation and critical exponents are extracted from numerical results in section Hv} section Ivl contains conclusions 
and perspectives. 



We consider the class of trivalent planar graphs, i.e all graphs without boundaries that can be drawn on a sphere 
and where each vertex is connected to exactly three neighbors. Moreover, two distinct vertices are connected by at 
most one link and no vertex is connected to itself. An example of <£> 3 graph is shown in figure [TJ 

For a given number of vertices N, all graphs can be obtained by gluing together N trivalent vertices in all possible 
ways satisfying the strong constraint of planarity. The Euler characteristic \ — N — Ni + Nf, where N, Ni, and Nf 
are, respectively, the number of vertices, links and faces, is fixed by the topology so that all planar graphs with the 
topology of a sphere share the same value x — 1. One more constraint, arising from the fact that each link is bounded 
by two vertices and that three links intersect at each vertex, reads 2Ni — 3N. These constraints imply that Nf and 
Ni are fixed if N is fixed. 

By giving each graph G a Boltzmann weight w(G), this class is turned into a statistical model of random graphs. 
In this paper, we choose w(G) — 1 for each graph. The partition function is given by 



where the sum is over all $ 3 graphs with N vertices as defined above. C(G) is a symmetry factor which avoids the 
overcounting of some symmetric graphs and is almost always equal to one for large graphs. There is an apparent 
regularity in these graphs. Namely connectivity, number of vertices, links and faces are fixed so that these graphs 
locally look like the honeycomb lattice. Moreover, the average size of a face, i.e the number of links surrounding a 
face, is equal to six for large N. However, the size of a face is, in fact, a random variable following a distribution P(Q) 
exactly known EH (see figure GJ), and the correlations associated to the variable Q are long range, decreasing as 
1/r 2 where r is the geodesic distance @. Note that for a face (i) with size Qi, the deviation of Qi from the average 
value < Q >— 6 defines the curvature R4 oc (Q, — 6)/Qi. 

This implies that these graphs are very different from the honeycomb lattice at large distance. For instance, there 
are objects called baby universes that induce a self-similar structure for these graphs. A baby universe is a connected 
part of a graph linked to it by a boundary called a neck. If we call B the size of a baby-universe and I the linear 
dimension of its neck, we must have I <C -B 1 / 2 , i.e. each baby universe is connected to the rest of the graph by a 
small boundary. On trivalent graphs, the minimum size of a neck is I = 3, and baby universes with such necks are 



II. RANDOM $ 3 PLANAR GRAPHS 
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Figure 2: Distribution of the size of the faces 

called minimum neck baby universes (minbus) . There are more of minbus than of the other baby universes and their 
distribution has been exactly calculated [12]. Moreover, baby universes can grow on other baby universes, so that 
each graph is like a tree made of a central part (the root) on which several branches of baby universes are growing. 
The global spherical topology ensures that branches do not intersect. This tree has a self-similar structure that can 
be seen as follows: consider the family of baby universes with a given ratio / = lj\fB between the linear size I of 
their neck and the square root of their size B. In a way, / defines a family of baby universes with the same shape. 
Then, let a{B — > 2B, f) be the (average) fraction of the total area of a graph included in baby universes with a shape 
characterized by / and with size ranging from B to 2B. It can be shown that a(B — » 2B, f) is independent of B. So, 
if we probe a graph, it looks similar at any scale, at least in the limit of large graphs. 

The fractal structure of these graphs described above is characterized by a Hausdorff (fractal) dimension du = 4 
This dimension is defined by the scaling law 

(N r ) ~ r dH 

where (N r ) is the average number of vertices at (geodesic) distance lower than r from any arbitrary vertex. 

III. TARGETED ATTACK AND ROBUSTNESS OF PLANAR $ 3 RANDOM GRAPHS 

A. Definitions 

A targeted attack is a process that systematically removes highly connected nodes of a graph. In fact, as each 
vertex of a $ 3 graph is always trivalent, we define targeted attack as nonrandom removal of faces - instead of vertices 
- with large number of edges. Note that, as <I> 3 planar graphs are dual to planar triangulations, this is equivalent to 
nonrandom removal of highly connected vertices of a triangulation. 

We introduce a cutoff number k max and define the removal process as follows: for a given graph G, we list all faces 
whose size - defined as the number of their edges - is greater than k max . Then, each face in this list is removed, i.e. all 
vertices and links which belong to its edges are removed. The resulting amputed graph G is generally not connected. 
Instead, G consists in a collection of connected smaller planar graphs called connected clusters. As we consider the 
statistical ensemble $ 3 |jv of planar $ 3 random graphs with N vertices, we are interested in average properties under 
targeted attacks, namely the properties of a graph G associated to a graph G picked at random in <& 3 |at. 

Robustness of a graph to targeted attacks means its ability to preserve its large scale connectivity after a nonrandom 
deletion of its highly connected nodes. This can be measured by its diameter, defined as the average geodesic 
distance between two arbitrary nodes, i.e. the average length of the shortest path between two arbitrary nodes 
Alternatively, as the effect of targeted attacks is to fragment graphs into smaller connected clusters, robustness can 
be measured by the average density of the largest connected cluster denoted P in this paper. Graphs in $> 3 \n are 
connected, so that P is equal to one before any targeted attack. To get more insight on the consequences of targeted 
attacks, we also measure the density of the second largest cluster, denoted S2, and the mean size of finite clusters, 
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Figure 3: Density of the largest cluster versus kmax- 
denoted S. The exact definition of S is the following [l4| 

s = s2 n ( s )/Yl s n ( s ) 

s s 

where n(s) is the average number of clusters made of s vertices and J2 S means that the largest cluster is excluded 
from the sum. 

B. Numerical experiment 

We generated graphs in Q 3 \n following the method described in appendix[Al For each graph in a Monte-Carlo series, 
we performed nonrandom removal of faces - for a given value of k max - and measured the size of the largest connected 
cluster by a method described in appendix[B] We simulated graphs in <£> 3 |at with N ranging from 1000 to 32000, and 
the parameter k max ranging from 8 to 30 for each value of AT, except for N = 32000 where k max was ranging from 
12 to 18. For each couple (N,k max ), we made 512 measurements separated by 1000 Ni local updating moves (flips) 
called Ti moves (see appendix [B]) . Errors were estimated using the standard jackknife method and error bars are 
systematically plotted in the figures. 



Results 



The density of the largest cluster P is plotted as a function of the cut-off k n 

P 1 rrr^T>l>o ™ ff.3 | 



large k max , P — 1 as graphs in $ d |jv are connected. We can see that for values of k 
k m ax, independently of the sizes of the graphs. Then, there is a fast fall of P dowi 
and more pronounced when the size is increased. This means that a targeted attack with k 



in figure [3l As expected, for very 
■nax > 17, P slowly decreases with 



greater than about 



17 does not significantly affect the connectivity of the graphs - i.e., for these values of k max , more than 70% of the 
vertices belong to the giant cluster. Conversely, k max has to be smaller than about 13 to completely fragment the 
graphs into small clusters. 

The average fraction of removed vertices is called x rm . It is a smooth decreasing function of k max as can be seen 
in figure IH It is interesting to plot P against x rm . This is shown in figure [5l 

We can see that a targeted attack significantly affects the connectivity only when a large fraction of the vertices 
are removed. For example, P > 0.7 for x rm ^ 0.13 and P > 0.5 for x rm < 0.17. In this sense, we can say that planar 
<I> 3 random graphs are rather robust to targeted attacks. 

The analysis of the density of the second largest cluster S2 and of the average size of finite clusters S, plotted 
in figures [6] and [7] respectively, gives more information on the mechanism of fragmentation of graphs under targeted 
attacks. For large k max , s 2 — 0, meaning that most of vertices belong to the largest cluster. Then, s 2 sharply increases 
when k max decreases and eventually reaches a peak for k max — 15. This is due to a small gradual decrease of the 
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Figure 4: Fraction of removed vertices versus k„ 
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Figure 5: Density of the largest cluster versus the fraction of removed vertices 



connectivity, i.e. a small fraction of vertices get disconnected from the largest cluster. However, this phenomenon is 
marginal as S2 does not exceed 0.15 and S is not greater than about 8% of the size. So, in this region k max > 15, the 
largest cluster completely dominates with only a few number of much smaller clusters. This is corroborated by the 
behavior of S in this region. As the largest cluster is excluded from this quantity, S ~ for very large k max . As k max 
decreases, S is essentially influenced by the second, third, etc, largest clusters and so it increases and reaches also a 
peak. Then, when k max continues to increase, i.e. when S2 and S go through their peak, there is a dramatic change 
in the fragmentation process. P undergoes a fast drop, meaning that there is an acceleration in the fragmentation 
and more and more vertices get disconnected from the initial largest cluster. The behaviors of S2 and S are similar. 
This is due to the fact that more and more small clusters appear, so that their average size is small, even for the 
largest, second largest, third largest, etc, clusters. 



IV. NON-UNIFORM PERCOLATION 



A. Percolation threshold 

The fragmentation process described above is very similar to the percolation mechanism This is, in fact, ex- 
pected as nonrandom removal of nodes under targeted attack is equivalent to percolation with non-uniform occupation 
probability [lj. In the case studied here, a vertex is occupied if and only if it does not belong to an edge of a face 
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Figure 6: Density of the second largest cluster versus fe^ 




Figure 7: Average size of finite clusters versus fe maa 



whose size is greater than k max - The behavior of P, S2 and S can then be interpreted by saying that there exists a 
critical value k c ~ 15 so that percolation takes place for k > k c . In the following, we consider k max as a continuous 
parameter; more precisely, it is as if the observables P, S2 and S were defined for all real (positive) values of k max 
but that numerical data were available only for integer values of k max . Note that k max can be seen as a cutoff on 
the local scalar curvature. By taking an appropriate continuum limit when N — > oo , the curvature is, of course, a 
continuous observable so that a cutoff on curvature has to be continuous in this limit. In other words, in any process 
of coarse graining, k max would effectively become a continuous parameter. 

In order to get a more precise value of k c , we define a finite size critical parameter k c (N) as the value of k max for 
which S2 is maximum for N fixed. By fitting the data of S2 with a Gaussian curve, we obtained k c (N). The result 
is shown in figure [8l The behavior of k c (N) can be estimated by standard finite size scaling analysis (see appendix 
ICj) . This allows us to predict the behavior of k c (N) approaching k c as k c (N) — k c ~ jq- 1 / l/d H w ^ n v > q a constant 
explained in section ITVB1 So, we fitted the values of k c (N) with the law k c (N) = k c + c N^ 1 / 1 " 1 " , with c a constant. 
We obtained k c = 15.8(2). 



B. Critical exponents z/, f3 and 7 



The point of view of non-uniform percolation with a percolation threshold k c strongly suggests that the observables 
S and S2 should obey the following scaling laws for k max near k c (13 . [H| (recall that we consider here k max as a 
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Figure 8: Finite size critical parameter k c (N) 
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with F and G two scaling functions. The exponent v characterizes the correlation length associated with a perco- 
lation transition. In particular, this exponent controls the way k c (N) tends to its infinite size value k c . So, the fit 
performed in the previous section also allows us to extract the value of v. The result is 1/i/djj = 0.4(2). In order to 
extract the values of exponents and 7, we have measured the maximum of S2 and S , which should scale as 



s™ ax {N) 
S max {N) 

The results are shown in figures l9l and flOl We fitted the data and obtained (3/vdH = 0.105(5) and 7/Wff = 0.77(1) 
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Figure 9: Maximum of density of the second largest cluster versus N, in logscale. The dotted line is the best fit 
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Figure 10: Maximum of the average cluster size versus N. The dotted line is the best fit 



C. Systematic errors 

Because of the discrete nature of the parameter k max , there are no available data for non- integer values of this 
parameter. This inevitably induces systematic errors - in particular, values of k max for which observables are maximum 
- certainly underestimated in the standard analysis used here. 



D. Comparing with uniform percolation exponents 

Critical exponent are known exactly for uniform percolation on planar <I> 3 random graphs fijvdn — 0.125, 
j/vcIh = 0.75 and l/vdn = 0.25. It is tempting to compare them with the values found in this work for non-uniform 
percolation fl/vdn = 0.105(5), "f/vdn = 0.77(1) and l/vda — 0.4(2). In view of systematic errors underestimated, 
it appears that values of exponents (3 and 7 are roughly compatible for uniform and non-uniform percolation. The 
exponent v seems very different from the exact value of uniform percolation. However, if we compare it with the 
value l/vdu = 0.489(9) found in numerical simulations Q, it is once again roughly compatible. In fact, as in uniform 
percolation, the fractal dimension dn is very sensitive to finite size effects. However, measures done in this work are 
not precise enough to conclude that critical exponents are the same for uniform and non-uniform percolation for this 
problem. It is not excluded either and this would imply a kind of universality between uniform and non-uniform 
percolation on these graphs. 



V. CONCLUSIONS 



In the first part of this work, the behavior of planar $ 3 random graphs under targeted attacks is investigated. The 
method consists in removing all vertices that belong to an edge of a face whose size is greater than a parameter k max . 
It appears that this family of graphs is rather robust to this nonrandom removal as k max has to be rather small 
(k m ax — 15 ) to significantly affect the global connectivity of these graphs, measured by the density of the largest 
cluster. This threshold corresponds to the removing of about 20% of the vertices. Note that the case of random 
failure is equivalent to uniform site percolation. This problem has already been investigated and the high value of the 
percolation threshold means that these graphs are rather fragile against random failure of nodes. 

In a second part, we consider the targeted attacks problem as non-uniform percolation. By measuring two ob- 
servables characterizing percolation clusters, namely the density of the second largest cluster S2 and the mean size of 
connected clusters S, and by extrapolating values of k max to non-integer values, we extract a critical value k c = 15.8(2) 
that can be interpreted as a percolation threshold for targeted attacks. It is also possible to look for scaling laws 
of observables in the vicinity of k max , by analogy to the uniform percolation problem. This allows us to measure 
(effective) standard critical exponent (3, 7 and v. Taking into account underestimated systematic errors due to the 
fact that k max is, in fact, a discrete parameter, critical exponents found here appear roughly compatible with their 
counterparts of uniform percolation. This is a very interesting fact as it would show that non-uniform and uniform 
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percolations are in the same class of universality. If true, this property could be a result of the hierarchical nature of 
these graphs that makes them look like trees of baby universes (BUs). 

Indeed, consider a vertex in a BU. For uniform percolation, a necessary condition for this vertex to belong to a 
giant cluster is that at least one vertex in the neck of the BU belongs to this giant cluster. The probability for this 
condition is rather small as necks are very small regions of the graph, so that the probability for a vertex to belong 
to a giant cluster is smaller if the vertex is in a BU. Consider now the case of targeted attacks. The probability for 
a vertex to be removed is, in fact, greater if it is in a BU as, on average, there is more curvature in BUs. Then, the 
probability that a vertex is in the boundary of a large face is greater if the vertex is in a BU. So, in both cases (uniform 
percolation and targeted attacks), the structure of BUs induces the same effect: the probability to be occupied is 
smaller for a vertex in a BU. However, this conjecture of universality between uniform and non-uniform percolation 
has to be more extensively and precisely investigated. 

Appendix A: GENERATING AND SAMPLING $ 3 GRAPHS 

The Monte-Carlo sequence starts with a graph in $ 3 |at randomly generated as follows: we start from a tetrahedron 
and, then, add vertices one by one in randomly chosen faces (triangles). Each new vertex is linked to the vertices 
of the corresponding triangle. The process is repeated until we obtain a graph (polyhedron) with N triangles. This 
graph becomes a <I> 3 graph - denoted G„ - by duality, i.e. by replacing each triangle of the polyhedron by a vertex 
linked to the vertices replacing the adjacent faces of the initial triangle. The resulting graph is topologically equivalent 
to a sphere with N trivalent vertices. Then, a series of graphs is obtained by u sing standard flips of links (T\ moves 
shown in figure [TT]) performed on randomly chosen links. It has been shown [lOj that T\ moves define an ergodic 
transformation in the ensemble <I> 3 |iv , so that whatever the starting graph G„, we obtain a nonbiased sampling with 
this method. 




Figure 11: Ti move involving faces 1,2,3 and 4. 



Appendix B: CLUSTER STRUCTURE 

In order to measure the connected cluster structure of graphs after a targeted attack, we use a breadth-first search 
algorithm similar to the Wolff algorithm [la ]. It recursively constructs all clusters for a given graph as follows: at 
"time" n, clusters ci,C2, . . . , c„ have already been detected. All sites in these clusters are labeled "visited". At time 
n+ 1, a not yet visited site (called v ) is chosen. By definition, it does not belong to any cluster already detected. v Q 
is the root of the new cluster c„+i; it is then labeled "visited" and put in a (empty) list Q. The following procedure 
is now recursively applied to Q : for each site v in Q, all not yet visited neighbors of v (on the graph) are added to 
Cn+i, labeled "visited" and put in Q whereas v is removed from Q . The procedure stops when Q is empty. Then, 
cluster c n +i is completely constructed. The algorithm stops when all sites have been visited. 

Appendix C: FINITE SIZE SCALING 

Finite size scaling analysis supposes that finite size corrections to scaling laws near a critical parameter k c are 
encoded by scaling functions depending on the ratio between the linear size L of the system and the correlation length 
£ . When I »(, finite size effects should not affect the system and scaling laws undergo no correction. On the 
contrary, when £ and L become comparable, finite size corrections of scaling laws are expected. For an observable O 
depending on a parameter k, an infinite size scaling law O ~ (k — k c ) z can be rewritten O ~ £, z ' v as the correlation 
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length scales as £ ~ (k — k c ) " ' . In particular, £ ~ L implies O ~ L 2 /". This can be summarized by the following law: 
O ~ L z l v F((k - fc c ) L 1 /") with F(x) a (scaling) function of the dimensionless ratio x = (k - k c )L x l v ~ (L/£,) 1/u . 
In order to interpolate between the cases L ^> £ and L ~ £, F(x) must verify F(x) — > 1 for x ~ 1 and F(x) — > x~ 2 for 
x » 1. Finite size scaling is a powerfull tool for extracting critical exponents by studying the variations of observables 
with the size of the system. This also provides a natural finite size critical parameter k c (L) : suppose that F(x) 
reaches a maximum for x = x Q . Then, for fixed L, the value of k giving a maximum for O obeys the following relation 
(k — k c ) L x l v — x a or, equivalently, k c (L) approaches k c when the size becomes infinite as k c (L) — k c ~ L~ x l v [l4l |. 
However, graphs in <& 3 |at have no explicit linear size, but the quantity N 1 ^" , where du is the Hausdorff dimension, 
plays this role so that finite size scaling laws are O ~ N z / udH F ((fc — k c ) N x / vdH ) and k c (N) approaches k c as 
k c (N) - k c ~ N- 1 ' 1 " 1 ". 
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